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1. Introduction 

In the last few years, lattice QCD simulations have made a substantial progress in controlling 
the systematic effects present in the determination of several important physical quantities allowing 
therefore for their direct contact with experiment (see ^ for recent reviews). Simulations in- 
cluding the light-quark flavours in the sea, as well as the strange and recently also the charm, with 
pseudoscalar masses below 300 MeV, lattice extents L > 2.0 fm and lattice spacings smaller than 
0.1 fm are currently being performed by several lattice groups. Such simulations will eventually 
allow for an extrapolation of the lattice data to the physical point and to the continuum limit while 
keeping also the finite volume effects under control. 

The European Twisted Mass collaboration (ETMC) has performed large scale simulations with 
Nf = 2 flavours of mass degenerate quarks using Wilson twisted mass fermions at maximal twist. 
Four values of the lattice spacing ranging from 0.1 fm down to 0.051 fm, pseudoscalar masses 
between 280 and 650 MeV as well as several lattice sizes (2.0 — 2.5 fm) are used to address the 
systematic effects. 

The light pseudoscalar meson is an appropriate hadron for investigating the systematic effects 
arising from continuum, infinite volume and chiral extrapolations, because its mass and decay 
constant can be obtained with high statistical precision in lattice simulations. Moreover, chiral 
perturbation theory (^PT) is best understood for those two quantities. As a consequence of this 
study one can extract other quantities of phenomenological interest, such as the u,d quark masses, 
the chiral condensate or the low energy constants of ;tPT. First results for the pseudoscalar mass 
mps and decay constant /ps from these Nf = 2 simulations can be found in Refs. [^-^. 

ETMC is currently performing A/^f = 2 + 1 + 1 simulations including in the sea, in addition 
to the mass degenerate light u,d quark flavours, also the heavier strange and charm degrees of 
freedom. Some first results for the pseudoscalar mass and decay constant from this novel setup 
were presented in p^ ]. 

In the following we will concentrate on the analysis of the Nf = 2 data for mps and /ps. 



2. Lattice Action and Setup 

In the gauge sector we employ the tree-level Symanzik improved gauge action (tlSym) [|]J 
The fermionic action for two flavours of maximally twisted, mass degenerate quarks in the so-called 
twisted basis [|l^, [l3|] reads 



Stm = [D[U]+ mo + 1^,75^'] Z(x)} , (2.1) 

where mo is the untwisted bare quark mass tuned to its critical value merit, i^q is the bare twisted 
quark mass, is the third Pauli matrix acting in flavour space and D\U] is the Wilson-Dirac 
operator. 

At maximal twist, i.e. mo = menu physical observables are automatically 0(a) improved with- 
out the need to determine any action or operator specific improvement coefficients (for a 



review see Ref. |14]). With this being the main advantage, one drawback of maximally twisted 
mass fermions is that flavour symmetry is broken explicitly at finite value of the lattice spacing, 
which amounts to O(a^) effects in physical observables. 
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Ensemble 


i8 


a [fm] 




mpsL 




mps [MeV] 


Di 


4.20 


0.051 


48^-96 


3.6 


0.0020 


280 


D2 






32^ • 64 


4.2 


0.0065 


510 


Ci 


4.05 


0.063 


32^ • 64 


3.3 


0.0030 


320 


C2 








4.6 


0.0060 


450 


C3 








5.3 


0.0080 


520 


C4 








6.5 


0.0120 


630 


C5 






243-48 


3.5 


0.0060 


450 


Ce 






203-48 


3.0 


0.0060 


450 


Bi 


3.90 


0.079 


243-48 


3.3 


0.0040 


330 


B2 








4.0 


0.0064 


420 


B3 








4.7 


0.0085 


480 


B4 








5.0 


0.0100 


520 


B5 








6.2 


0.0150 


640 


Be 






323 . 64 


4.3 


0.0040 


330 


Bi 








3.7 


0.0030 


290 


A2 


3.80 


0.100 


243-48 


5.0 


0.0080 


410 


A3 








5.8 


0.0110 


480 


A4 








7.1 


0.0165 


580 



Table 1: Ensembles with N{ = 2 dynamical flavours produced by die ETM collaboration. We give the 
ensemble name, the values of the inverse bare coupUng j3 = 6/^0, an approximate value of the lattice spacing 
a, the lattice volume V = T in lattice units, the approximate value of nip^L, the bare quark mass in 
lattice units and an approximate value of the light pseudoscalar mass mps. 



For details on the setup, tuning to maximal twist and the analysis methods we refer to Refs. [||, 
^ |6|]. Recent results for light quark masses, meson decay constants, the pion form factor, n-n 
scattering, the light baryon spectrum, the rj' meson and the co — p mesons mass difference are 
available in Refs. [|l|, |l|, ?, [Hl-ii- 

Flavour breaking effects have been investigated for several quantities ||3|, ^ ^ ^, |T7|]. With the 
exception of the splitting between the charged and neutral pion masses, other possible splittings 
so far investigated are compatible with zero. These results are in agreement with a theoretical 
investigation using the Symanzik effective Lagrangian [ pT} , p2[ ]. 

A list of the Nf = 2 ensembles produced by ETMC can be found in table [T[ 



3. Results 

3.1 Scaling to the Continuum Limit 

Here we analyse the scaling to the continuum limit of the pseudoscalar meson decay constant 
fps at fixed reference values of the pseudoscalar meson mass mps and of the lattice size L (we refer 
to [Q, ^, ^] for details). The aim of this scaling test is to verify that discretisation effects are indeed 
of 0{a^) as expected for twisted mass fermions at maximal twist. 
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Figure 1: (a) rg/ps as a function of (romps) (b) Continuum limit scaling: ro/ps as a function of {a/ro) at 
two fixed values of romps. 



In order to compare results at different values of the lattice spacing it is convenient to measure 



the hadronic scale ro/a [23]. It is defined via the force between static quarks at intermediate 
distance and can be measured to high accuracy in lattice QCD simulations. For details on how we 
measure ro/a we refer to Ref. [H]. 



In figure |l(a)| we plot the results for ro/ps as a function of (romps) . The vicinity of points 



coming from different lattice spacings along a common curve is an evidence that lattice artifacts 



are small for these quantities. This is indeed confirmed in figure |l(b)| where the continuum scaling 
of ro/ps is illustrated: the very mild slope of the lattice data shows that the expected 0{a^) scaling 
violations are small. The result of a linear extrapolation in [a/roY to the continuum limit is also 
shown. 

3.2 xVT Description of Finite Size Effects 

At the level of statistical accuracy we have achieved, finite size effects (FSE) for /ps and mps 
cannot be neglected. It is therefore of importance to study whether FSE can be described within 
the framework of chiral perturbation theory. This requires to compare simulations with different 
lattice volumes while all other parameters are kept fixed, like for instance ensembles C2, C5 and Cg 
or B\ and in table |l]. For all these ensembles mpsL > 3 holds, which is believed to be needed 
for ;^PT formulae to apply. Given the smallness of the lattice artifacts in /ps and nips (as discussed 
in the previous section), we proceed to compare the measured finite size effects to predictions of 



continuum %VY at NLO [ |24| ] (denoted GL) and in the form of the resummed Liischer formula as 
described in Ref. ^ (for short CDH). 

Wenote/?o = [0(L = oo)-0(L)]/0(L = oo)the relative FSEfor the observable 6>G {mps,/ps}. 
The results for R^^^\ R^^ and R^^^ are compiled in table |2[ We observe that the CDH formulae 
tend to provide an appropriate description of the lattice data. A more detailed description of FSE 
in our /ps and m^s data was presented in Ref. 
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a [fm] 


mpsLi - 


-> mpsL2 




R^^ [%] 


7?'^^^ [%] 










[Li ^L2] 


[Li ^ L2] 


[Li ^ L2] 


mps 


0.079 


3.3- 


-^4.3 


-1.8 


-0.4 


-1.0 


/ps 


0.079 


3.3- 


-^4.3 


+2.6 


+2.1 


+2.3 


mps 


0.063 


3.0- 


-^4.6 


-6.1 


-1.7 


-5.9 


/ps 


0.063 


3.0- 


-^4.6 


+10.7 


+6.3 


+8.5 



Table 2: Comparison of measured relative FSE, Rq, to estimates from ;ifPT formulae. 



3.3 xVT Description of the light-quark Mass Dependence 

The chiral extrapolation of lattice data down to the physical point is currently one of the main 
sources of systematic uncertainties in the lattice results. The possibility to rely on an effective 
theory such as ;^PT to perform this extrapolation is therefore of great importance to quote accurate 
results from lattice simulations. On the other hand, while smaller quark masses are being simulated, 
the possibility to perform a quantitative test of the effective theory as well as to measure the low 
energy parameters of its Lagrangian becomes more and more realistic. 

We shall now present the results of a combined chiral, infinite volume and continuum ex- 
trapolation of mps and /ps for two values of the lattice spacing (corresponding to j3 = 3.90 and 
j8 = 4.05). We use ro/c? to relate data from the two lattice spacings and a non-perturbative deter- 



mination of the renormalisation factor Zp [ 25 ] in order to perform the fit in terms of renormalised 
quark masses. This analysis closely follows those presented in Refs. [Q, |[, ^, ^] to which we refer 
for more details. 

We perform combined fits to our data for /ps, mps, ro/a and Zp at the two values of j3 with the 
formulae: 



?"o/ps = ''0/0 



{rmvsf = lA [1 + ^ log (^) + + A„ps(«/ro)^] 



(3.1) 



with ^ = 2BoHR/{4nfof , Xn = 2BoHr , IXr = jX^/Zp, /(, = V2Fo. T,^° denote the continuum 
NNLO terms of the chiral expansion [26], which depend on Ai_4 and and kf, and K^^^{L) 



the finite size corrections [27]. Based on the form of the Symanzik expansion in the small quark 



mass region, we parametrise in eq. (3. 1 ) the leading cut-off effects by the two coefficients Dy-p^ ,„pg. 
Setting £)/ps,mps = is equivalent to perform a constant continuum extrapolation. Similarly, setting 

rNNI 



rNNL° = corresponds to fit to NLO Z^T- 



From the fit parameters coming from the quark mass dependence predicted by ;fPT (in par- 
ticular from A3 4, Bo and /o) the low energy constants £3 4 and the chiral condensate £ can be 
determined. 

By including or excluding data points for the heavier quark masses, it is in principle possible 
to explore the regime of masses in which NLO and/or NNLO SU(2) ;t^PT applies. We have actually 
generalised this procedure in order to estimate all the dominant sources of systematic uncertainties 
that can be addressed from our setup, which include, discretisation effects, the order at which we 
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Figure 2: Quark mass dependence: (a) Data for (romps)^ / rojiR as a function of ro/Xfi. (b) Data for ro/ps as a 
function of ro/ifl. The data are from ensembles B1.2, 3,4,6 and Ci 2.3.5 and the fit to this data is of type B. Note 
that in these figures we did not propagate the errors of ro and Zp. 



work in ;^PT or finite size effects. The idea is to use different fit ansatz (see below) on a given 
data-set and to repeat this same procedure over different data-sets: by weighting all these fits by 
their confidence level we construct their distribution and estimate the systematic error from the 
associated 68% confidence interval. 
The fit ansatze we consider are: 

• Fit A: NLO continuum x^T, r„^/^° = 0, D,„p^jp^ = 0, priors for roAi,2 

• Fit B: NLO continuum x^T, = 0, D,„p^jp^ fitted, priors for roAi,2 

• Fit C: NNLO continuum ;^PT, A)ips,/ps = 0' priors for tqAi 2 and ^m,f 

• Fit D: NNLO continuum ;^PT, £)mps./ps fitted, priors for tqAi 2 and kM.F 

The choice of the different data-sets (each of them including data for different lattice spacings, 
quark masses and physical volumes) is made in order to quantify how the quality of the fit is 
modified when including/excluding data from e.g., a given mass region or with a given volume. 
The data-sets considered in the fits are listed in Ref. [pi]. 



As an example we show in figures [2(a)| and [2(b)| the result for a fit of type B on data-set 
composed of ensembles Bi,2,3,4,6 and Ci.2,3,5 

The main physical results we obtain from this analysis are the light quark mass mJ^J(/x = 
2GeV) = 3.54(26)MeV, the pseudo scalar decay constant in the chiral hmit /o = 122(l)MeV, the 
scalar condensate [lF^{iJ. = 2GeV)]^/^ = 270(7)MeV and fj,/fo = L0755(94). We furthermore 
extract accurate values for other low energy constants of ;tPT, in particular £3 = 3.50(31) and 
£4 = 4.66(33). The errors are statistical and systematic errors summed in quadrature. 
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3.4 Discussion and Conclusion 

Here we collect a short list of observations coming from a set of ;^PT fits. A complete descrip- 
tion of a large set of combined fits, including the details on the estimates of the systematic effects, 
was presented in Ref. 

We observe that including in the fits pseudoscalar masses mps > 520 MeV decreases signifi- 
cantly the quality of the NLO fits (x^/dof » 1). This indicates that the applicability of NLO ;^PT 
in that regime of masses is disfavoured. 

On the contrary, extending the fit-range to a value of nips ~ 280 MeV preserves the good 
quality of the fit and gives compatible values for the fit parameters. This result makes us confident 
that the extrapolation to the physical point is trustworthy. 

Including lattice artifacts in the fits gives results which are compatible to those where Omps,/ps 
is set to zero but with a somehow better ;|^^/dof. We observe that the values of the fitting parameters 
A«ps,/ps compatible with zero within two standard deviations. This is in line with the small 
discretisation effects observed in the scaling test. 

The inclusion of NNLO terms produces similar results to the NLO fits in the quark mass 
region corresponding to mps E [280,520] MeV. When fitting data only in this mass region {i.e. 
when excluding from the fit the heavier masses at mps ~ 650 MeV), we observe that the fit curve 
at NLO lies closer to those data points (heavier masses) than the NNLO one. On the other hand, 
when including the heavier masses, the NNLO fit is able describe these data points but the quality 
of the fit is somehow reduced with respect to the NLO fit. To improve the sensitivity of our lattice 
data to ;^PT at NNLO, additional data points would be needed. 

We have presented determinations of /ps and nips and their continuum, infinite volume and 
chiral extrapolations. A complete description of the results was presented in Ref. [j8|]. 

We thank all members of ETMC for the most enjoyable collaboration. 
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